Uncertainties in nuclear transition matrix elements for neutrinoless /3/3 decay within 

the PHFB model 
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The nuclear transition matrix elements M'-"''-' for the neutrinoless double beta decay of ^'''^^Zr, 
^*'^°°Mo, lo^Ru, ^^°Pd, i28,i30r|.g isojv^j isotopes in the case of 0+ 0+ transition are calculated 
using the PHFB wave functions, which are eigenvectors of four different parameterizations of a 
Hamiltonian with pairing plus multipolar effective two-body interaction. Employing two (three) 
different parameterizations of Jastrow-type short range correlations, a set of eight (twelve) different 
nuclear transition matrix elements M^""^ is built for each decay, whose averages in conjunction with 
their standard deviations provide an estimate of the model uncertainties. 

PACS numbers: 21.60.Jz, 23.20.-g, 23.40.Hc 



I. INTRODUCTION 

Ascertaining the mass and nature of neutrinos requires 
the analysis of observational data obtained from three 
complementary experiments, namely single- /3 decay, neu- 
trino oscillation and neutrinoless double beta (/?/?) de- 
cay. In any gauge theoretical model with spontaneous 
symmetry breaking, the observation of (/3/?)oi, decay im- 
plies non-zero mass of Majorana neutrinos at the weak 
scale independent of the underlying mechanisms [T], 
The varied scope and far reaching nature of experimen- 
tal as well as theoretical studies devoted to (/3/3)o!y decay 
over the past decades have been excellently reviewed by 
Avignone et al. 31 and references there in. 

The observed experimental limits on the half-life 
of decay have already provided stringent limits on 

the associated gauge theoretical parameters [1] . The reli- 
ability of extracted gauge theoretical parameters depends 
on the accuracy of nuclear transition matrix elements 
(NTMEs). For a given transition, different NTMEs are 
obtained employing distinct nuclear models, and for a 
given model, they also depend on the model space and 
effective two-body interaction selected. Other uncertain- 
ties are related with the inclusion of pseudoscalar and 
weak magnetism terms in the Fermi, Gamow- Teller and 
tensorial NTMEs Q , finite size as well as short range 
correlations [7l-[l0|. and the use of two effective values of 
the axial- vector coupling constant qa- 

The spread between the calculated NTMEs provides a 
measure of the theoretical uncertainty [ll|. In the case 
of the well studied ^^Ge isotope, it was observed that 
the calculated decay rates differ by a factor of 6-7. The 
effective neutrino mass {m^) is inversely proportional to 
the square root of T^'j,^. Hence, the uncertainty in the 
effective neutrino mass is about 2 to 3. For example, 
from the experimental limit Tjyj > 1-6 x 10^^ yr 
the upper limits on (m^) range between 0.4 eV and 1.0 
eV, depending on the NTME [l3|-[l|. If the {PP)^^ de- 



cay were observed in several nuclei, the comparison of 
calculated ratios of the corresponding NTMEs-squared 
and the ratios of half-lives could also test the validity 
of nuclear structure calculations in a model independent 
way p^ . 

Rodin et al. [l^l have estimated the theoretical uncer- 
tainty employing two models, the QRPA and RQRPA, 
with three sets of basis states and three realistic two-body 
effective interactions. Different strategies to remove the 
sensitivity of QRPA calculations on the model param- 
eters have been proposed [H, [l^. Further studies on 
uncertainties in NTMEs due to short range correlations 
using the unitary correlation operator method (UCOM) 
and self-consistent coupled cluster method (CCM) Q 
have been also carried out by Faessler and coworkers. 

Up to now, the QRPA model and its extensions 
have been the most successful models in correlating the 
single-/? GT strengths and half-lives of P^)2v de- 
cay and the first in explaining the observed suppres- 
sion of A/21/ [lOillH- Nonetheless, the large scale shell 
model (LSSM) calculations of Strasbourg-Madrid group 
are quite promising [23|. Deformation has been included 
at various levels of approximation in the QRPA for- 
malism 23-25]. Recently, the effects of pairing and 
quadrupolar correlations on the NTMEs of (/3~/3~)qj, de- 
cay have also been studied in the interacting shell model 
(ISM) 1,11^ and the projected-Hartree-Fock-Bogoliubov 
(PHFB) model [13, [2|. 

The PHFB model, in conjunction with pairing plus 
quadrupole-quadrupole (PQQ) [1^ interaction has been 
successful in the study of the O"*" — ?> 0+ transition of 
{P~ P~)2u decay, where it was possible to describe the 
lowest excited states of the parent and daughter nu- 
clei along with their electromagnetic transition strengths, 
as well as to reproduce their measured /3/3 decay rates 
[30I [3l| . The PHFB model is unique in allowing the de- 
scription of the /3/3 decay in medium and heavy mass 
nuclei by projecting a set of states with good angular 
momentum, while treating the pairing and deformation 
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degrees of freedom simultaneously and on equal footing. 
On the other hand, in the present version of the PHFB 
model, the structure of the intermediate odd Z-odd N 
nuclei and hence, the single l3 decay rates and the distri- 
bution of GT strength can not be studied. Notwithstand- 
ing this limitation, it is a convenient choice to examine 
the explicit role of deformation on the NTMEs. In the 
study of /3~/3~ decay, there are four noteworthy qbser- 
vations in connection with deformation effects p7l . [28j , 
namely: 

(i) There exists an inverse correlation between the 
quadrupole deformation and the size of NTMEs 
M2., Af(o-) and M^j^^l 

(ii) The NTMEs are usually large in the absence of 
quadrupolar correlations; they are almost constant 
for small admixture of the QQ interaction and sub- 
stantially suppressed in deformed nuclei. 

(iii) In agreement with the observations made by 
Simkovic et al. [H, the NTMEs have a well de- 
fined maximum when the deformation of parent and 
daughter nuclei are similar, and they are quite sup- 
pressed when the difference in the deformation is 
large. 

(iv) The deformation effects are of equal importance in 
case of {P~ P~)2i, and (/3~/3~)oi/ decay. 

In earlier works, we have calculated NTMEs M21, for 
the {§_l3-)2^ [Him and M^^"'^ for the {P_J~)q^ de- 
cay [23] with the PQQ effective interaction |29l| , and the 
effect of hexadecapolar correlations (HH) [2^ on the cal- 
culated spectroscopic properties and (/3~/3~) decay rates 
has been studied. In the present work, we employ two 
different parameterizations of the QQ interaction, with 
and without the HH correlations. Further, the NTMEs 
jy^(Oi/) g^j.g calculated with three different parametriza- 
tions of Jastrow type of SRC employing the four sets of 
wave functions. The twelve NTMEs provide a reason- 
able sample for estimating the associated uncertainties. 
In Sec II, the PHFB formalism employed to describe the 
(/3~/3~)qj^ decay with the inclusion of the finite size of 
the nucleons and short range correlations is shortly re- 
viewed. In Sec. Ill, the four different parameterizations 
of the pairing plus multipole Hamiltonian are introduced, 
the calculated NTMEs vis-a-vis their radial evolution are 
analyzed, and their average values as well as standard 
deviations are estimated. Subsequently, the latter are 
employed to obtain upper limits on the effective mass of 
light Majorana neutrinos. Conclusions are given in Sec. 
IV. 



light neutrinos for the 0+ 0+ transition is given by 



Ti%(0+ ^ 0+) 



where the NTMEs mI""^ are given by 



(1) 
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with 



Of = 



— ) H{n2), OGT^<Ji-a2H{n2) (3) 
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and 



Hiru) 



R<i){Ari2) 



ri2 
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The origin of the neutrino potential H{ri2) is due to the 
exchange of light Majorana neutrinos between nucleons 
being considered as point particles. To take the finite 
size of nucleons into account, neutrino potential H(ri2) 
is folded with a dipole form factor and rewritten as 



Hiru 



where 



AttR 



3 exp(jq-ri2) 



A2 



N 



q{q + A) VA2 + g- 



{Ei + Ef) 



(5) 



(6) 



and the cutoff momentum A= 850 MeV '2?]. 

The short range correlations (SRC) are produced 
by the repulsive nucleon-nucleon potential generated 
through the exchange of p and lo mesons. They have been 
included in the calculations of Af'"") for the {P~ de- 
cay through the phenomenological Jastrow type of corre- 
lations with Miller-Spenser parametrization 35] , effective 
operators [36|], exchange of w-meson ^STi], UCOM [7, 3^ 
and self-consistent CCM Q. It has been observed that 
the effects due to the Jastrow type of correlations with 
Miller-Spenser parametrization are usually strong [36| . 
where as the UCOM and self-consistent CCM have weak 
effects. Further, Simkovic et al. Q have shown that it is 
possible to parametrize the SRC effects of Argonne V18 
and CD-Bonn two nucleon potentials by the Jastrow type 
of correlations within a few percent accuracy. Explicitly, 
the effects due to the SRC can be incorporated in the 
calculation of M^'^^'> through the prescription 



(7) 



II. THEORETICAL FORMALISM 

In the Majorana neutrino mass mechanism, the inverse 
half- life of the (/3^/?^)o^ decay due to the exchange of 



with 



fir) 



1 — ce 



(8) 



where a = 1.1, 1.59 and 1.52 /jti"^ b = 0.68, 1.45 and 
1.88 fm-^ and c ^ 1.0, 0.92 and 0.46 for Miller-Spencer, 
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Argonne V18 and CD-Bonn NN potentials, respectively. 
In the next section, the NTMEs M^°'^^ are calculated in 
the PHFB model by employing these three sets of param- 
eters for the SRC, denoted as SRCl, SRC2 and SRC3, 
respectively. 

The three functions f{r) are plotted in Fig. [TJ They 
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FIG. 1: Radial dependence of /(r) for the three different pa- 
rameterizations of the SRC. 

have similar forms, but differ in its value at the origin, 
and at the position of its maximum, which lies at 1.54, 
1.15 and 1.09 fm for SRCl, SRC2 and SRC3, respectively. 
They have clear influence on the radial evolution of the 
(/3~/3~)qj^ decay matrix elements discussed below. 

The calculation of M^^"^) in the PHFB model has been 
discussed in Ref. [l^ and one obtains the following ex- 
pression for NTMEs Aff ""^ of W~p-)o^ decay 



M. 



X / »^(z.^'),(z+2.^r-2)(^') ^ (a/3|Ofc|7(5) 

"/57<5 



E 



(A- 



Z+2,N-2 



1 + F^"' le^f^"'* 

^ ' ^ ZM\" IJ Z+2,N-2 



sinOdO, (9) 



where 



1/2 



det (^H-i^('^)/W^' 

1 /2 

det {l + F^'^V^''^^) d-l^^{e)sin{e)deiiQ) 



and 

''^{Z,N),(Z+2,N-2) (0) 



A^i ( 1 -X- W^^^l f ("^^ 

aei I 1 -t- ^ z,nJ Z+2,N-2 



aei. [ 1 -f- ^ z,nJ Z+2,N-2 



1/2 



1/2 
(11) 



The 7r(;^) represents the proton (neutron) of nuclei in- 
volved in the (/3~/3~)oj, decay process. The matrices fz,N 
and Fz,n{0) are given by 

i 



The extra factor 1/4 in the Eq. (28) of Ref. [23] should 
not be there. 



III. RESULTS AND DISCUSSIONS 

The model space, single particle energies (SPE's), pa- 
rameters of the PQQ type of effective two-body inter- 
action and the method to fix them have been already 
given in Refs. [27l[30ll3ll. Presentlv. we use the effective 
Hamiltonian written as [28| 

H = Hsp + V{P) + V{QQ) + V{HH), (14) 

where H^p, V{P), V{QQ) and V{HH) denote the 
single particle Hamiltonian, the pairing, quadrupole- 
quadrupole and hexadecapole-hexadecapole parts of 
the effective two-body interaction, respectively. The 
quadrupole-quadrupole part of the effective two-body in- 
teraction V{QQ) has three terms, namely the proton- 
proton, the neutron-neutron and the proton-neutron 

denoted by \ 

[13, S Eli, 

of the like particle components of the QQ interaction 
were taken as X2pp — X2nn — 0.0105 MeV h-^, where 
h is the oscillator parameter. The strength of proton- 
neutron component of the QQ interaction X2pn was var- 
ied so as to fit the experimental excitation energy of the 
2"*" state, E2+ ■ In the present work, we also employ an 
alternative isoscalar parametrization of the quadrupole- 
quadrupole interaction, by taking X2pp = X2nn = X2pn/2. 
In this case, the three parameters are varied together to 
fit i?2+. We will refer to these two parameterizations 
of the quadrupole-quadrupole interaction as PQQl and 
PQQ2. 

Employing either method, the experimental excita- 
tion energies of 2+ state E2+ [S^ can be reproduced 
within about 2% accuracy. The maximum change in 
i?4+ and Ek + energies with respect to PQQl interac- 
tion [sO, [3l[ is about 5% and 18%, respectively. The 
reduced B{E2:0^ — > 2+) transition probabilities, defor- 
mation parameters /?2, static quadrupole moments Q(2+) 



ones, whose coefficients are denoted by X2pp,X2rm and 
X2pn, respectively. In Refs. [13, HO, HH, the strengths 
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and gyromagnetic factors 5(2''") are in an overall agree- 
ment with the experimental data [iO, |4l| for both the 
parametrizations. In the case of PQQ2 parametriza- 
tion, the maximum change in the calculated NTMEs 
M21, for the 0+ — > 0+ transition with respect to PQQl 
parametrization is about 21% but for ^^Zr isotope. 

The HH part of the effective interaction V{HH) is 
given as [ 28j 



V{HH) 



x(/3|g4_^|(5) aj^aj, a^,, (15) 



with = r*y4^(0, 0). The relative magnitudes of the 
parameters of the HH part of the two body interaction 
are calculated from a relation suggested by Bohr and 
Mottelson [42] . The approximate magnitude of these con- 
stants for isospin T = is given by 



47r 



XA 



2A-H A(r2^-2) 



/or A = 1,2,3,4. 



(16) 



and the parameters for the T — \ case are approximately 
half of their T = counterparts. Presently, the value of 
X4 = 0.2442 X2^"^/^^"'' for T = 1, which is exactly half 
of the T = case. 

We refer to the calculations which include the hex- 
adecapolar term HH as PQQHH. We end up with 
four different parameterizations of the effective two- 
body interaction, namely PQQl, PQQHHl, PQQ2 and 
PQQHH2. 



instead of A, is of the order of 10 %. It confirms that 
the dependence of NTMEs on average excitation energy 
A is small for the (/3~/3~)q,^ decay and the validity of the 
closure approximation is quite satisfactory. 

The variation in M^'^'^^ due to the different parameter- 
izations of the Hamiltonians (presented in the different 
rows) lies between 20-25%. It is noticed in general but 
for i28Te isotope that the NTMEs evaluated for both pa- 
rameterizations of the quadrupolar interaction are quite 
close. The inclusion of the hexadecapolar term tends to 
reduce them by amounts which strongly depend on the 
specific nuclei. 

The inclusion of SRC in the approximation of point nu- 
cleons (P-l-S) induces an extra quenching in the NTMEs 
M^O'') , which can be of the order of 18-23% for SRCl, to 
negligible for SRC3. The dipole form factor (F) always 
reduces the NTMEs by 12-15% in comparision to the 
point particle case. Adding SRC (F-f S) can further re- 
duce the transition matrix elements, for SRCl, or slightly 
enhance them, partially compensating the effect of the 
dipole form factor. It is interesting to note that the ef- 
fect of F-SRC2 is almost negligible, i.e., nearly the same 
as F. 

The radial evolution of M^'^") has been studied in the 
QRPA by Simkovic et al. 7] and in the ISM by Menendez 
et al. [43] by defining 



(17) 



A. SRC and radial evolutions of NTMEs 

In Table U the NTMEs M'-""^ evaluated using 
the HFB wave functions in conjunction with PQQl, 
PQQHHl, PQQ2, PQQHH2 interactions and three dif- 
ferent parametrizations of the Jastrow type of SRC for 
the nuclei ^^'^^Zi, 98,100^0, 104ru^ nop^^ 128,130^6 and 
-'^^"Nd are displayed. The average energy denominator 
A has been taken as ^ = 1.12^^/^ yi^y following Hax- 
ton's prescription [l^ . The NTMEs are calculated in the 
the approximations of point nucleons (P - 2nd and 3rd 
columns), finite size of nucleons (F - 7th column), point 
nucleons with SRC (P+S - 4th to 6th columns), and fi- 
nite size plus SRC (F+S - last three columns). To obtain 
additional information on the stability of the estimations 
of NTMEs M^"'"), they are also calculated for A/2 in the 
energy denominator in the case of point nucleons, given 
in the column 3. 

We present the relative changes in NTMEs Af (0") ( m 
%) due to the different approximations in Table HT] In 
each row, i.e. for each set of wave functions, the reference 
NTMEs M^"") are those calculated for point nucleons 
without SRC, given in the second column of Table ID 
It can be observed that the relative change in NTMEs 
M(o^), when the energy denominator is taken as A/2 



In both QRPA and ISM calculations, it has been estab- 
lished that the contributions of decaying pairs coupled to 
J — and J > almost cancel beyond r 3 fm and the 
magnitude of C'-^"-' for all nuclei undergoing (/3~/3~)gj^ 
decay are the maximum about the internucleon distance 
r « 1 fm. In Fig. [21 we plot the radial dependence 
of (7*^°"^ due to PQQl parametrization of the effective 
two body interaction for six nuclei, namely ^^Zr, ^™Mo, 
""Pd, i^saao^g j^jj^j 150^^^ The radial evolution of M^O"^) 
is studied for eight cases, namely P, P-l-SRCl, P-I-SRC2, 
P+SRC3, F, F-f SRCl, F-)-SRC2 and F-I-SRC3. In addi- 
tion, the effects due to the finite size and SRC are made 
more transparent in Fig. [3] by plotting them for different 
combinations of P, F and SRC. In case of point nucle- 
ons, it is noticed that the (7*^°"^ are peaked at r = 1.0 fm 
and with the addition of SRCl, the peak shifts to 1.25 
fm. However, the magnitude of C^-^"^^ are increased for 
SRC2 and SRC3 with unchanging peak position. In the 
case of ENS, the C*^'^"-' are peaked at r = 1.25 fm, which 
remains unchanged with the inclusion of SRCl, SRC2 
and SRC3. However, the magnitudes of (7*^°"^ change in 
the latter three cases. The above observations also re- 
main valid with the other three parametrizations of the 
effective two-body interaction. 
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TABLE I: Calculated NTMEs M^^"^ in the PHFB model with four different parameterizations of the effective two-body 
interaction, three different parameterizations of SRC, with nucleons taken as point particles (P) or with a dipole form factor 
(F), for the (/3-/3-)o^ decay of ''^■'■"'Zv, 9»>iooMo, ^o^u, "^Pd, i28>i30Te and ^^ONd isotopes. 

Nuclei 



P P+S F F+S 









A 


A/2 


SRCl 


SRC2 


SRC3 






SRCl 


SRC2 


SRC3 


y^Zr 


PQQl 


5 


4382 


5.8729 


4 


3021 


5 


0644 


5 


4097 


4 


6891 


4 


0690 


4 


6639 


4 


8383 




PQQHHl 


5 


0015 


5.3947 


3 


9472 


4 


6528 


4 


9734 


4 


3069 


3 


7315 


4 


2820 


4 


4441 




PQQ2 


5 


1183 


5.5700 


4 


1781 


4 


8201 


5 


1063 


4 


4912 


3 


9802 


4 


4818 


4 


6259 




PQQHH2 


4 


7626 


5.1356 


3 


7492 


4 


4266 


4 


7348 


4 


0955 


3 


5424 


4 


0708 


4 


2266 


^^Zr 


PQQl 


3 


9517 


4.2741 


3 


0829 


3 


6622 


3 


9257 


3 


3828 


2 


9068 


3 


3590 


3 


4923 




PQQHHl 


3 


9363 


4.2413 


3 


0330 


3 


6333 


3 


9072 


3 


3459 


2 


8507 


3 


3192 


3 


4578 




PQQ2 


3 


7804 


4.0875 


2 


9450 


3 


5017 


3 


7549 


3 


2335 


2 


7758 


3 


2103 


3 


3385 




PQQHH2 


3 


7035 


3.9891 


2 


8470 


3 


4158 


3 


6753 


3 


1442 


2 


6745 


3 


1182 


3 


2497 


''^Mo 


PQQl 


8 


7743 


9.5345 


7 


0859 


8 


2350 


8 


7484 


7 


6507 


6 


7322 


7 


6297 


7 


8884 




PQQHHl 


8 


1669 


8.8420 


6 


5385 


7 


6442 


8 


1395 


7 


0846 


6 


1984 


7 


0618 


7 


3114 




PQQ2 


8 


8254 


9.5866 


7 


1202 


8 


2806 


8 


7992 


7 


6907 


6 


7630 


7 


6695 


7 


9307 




PQQHH2 


8 


0911 


8.7589 


6 


4723 


7 


5712 


8 


0636 


7 


0154 


6 


1344 


6 


9925 


7 


2406 


iooMo 


PQQl 


8 


5939 


9.2939 


6 


8691 


8 


0522 


8 


5763 


7 


4413 


6 


5036 


7 


4282 


7 


6920 




PQQHHl 


8 


2130 


8.8577 


6 


5174 


7 


6763 


8 


1915 


7 


0822 


6 


1597 


7 


0654 


7 


3248 




PQQ2 


8 


6571 


9.3633 


6 


9212 


8 


1116 


8 


6391 


7 


4972 


6 


5534 


7 


4838 


7 


7493 




PQQHH2 


7 


4186 


7.9968 


5 


8774 


6 


9312 


7 


3994 


6 


3904 


5 


5520 


6 


3756 


6 


6113 


lo^'Ru 


PQQl 


6 


2757 


6.7734 


4 


9743 


5 


8753 


6 


2705 


5 


4007 


4 


6942 


5 


3989 


5 


5975 




PQQHHl 


5 


7976 


6.2339 


4 


5484 


5 


4102 


5 


7895 


4 


9596 


4 


2809 


4 


9548 


5 


1454 




PQQ2 


5 


9034 


6.3698 


4 


6777 


5 


5267 


5 


8989 


5 


0789 


4 


4137 


5 


0777 


5 


2647 




PQQHH2 


5 


3786 


5.7803 


4 


2143 


5 


0176 


5 


3711 


4 


5974 


3 


9648 


4 


5931 


4 


7708 


"Opd 


PQQl 


10.1361 


11.0441 


8 


1250 


9 


5068 


10.1167 


8 


7918 


7 


6982 


8 


7783 


9 


0850 




PQQHHl 


8 


5617 


9.2893 


6 


7742 


7 


9988 


8 


5408 


7 


3694 


6 


3963 


7 


3535 


7 


6262 




PQQ2 


9 


7208 


10.5944 


7 


7929 


9 


1163 


9 


7011 


8 


4328 


7 


3842 


8 


4187 


8 


7128 




PQQHH2 


9 


0246 


9.8138 


7 


1864 


8 


4447 


9 


0023 


7 


7985 


6 


7982 


7 


7816 


8 


0621 


128m 


PQQl 


4 


3415 


4.7394 


3 


4372 


4 


0474 


4 


3219 


3 


7417 


3 


2499 


3 


7258 


3 


8639 




PQQHHl 


4 


8152 


5.2111 


3 


7261 


4 


4626 


4 


7931 


4 


0916 


3 


4994 


4 


0740 


4 


2401 




PQQ2 


5 


1422 


5.6212 


4 


1056 


4 


8082 


5 


1233 


4 


4521 


3 


8893 


4 


4374 


4 


5956 




PQQHH2 


5 


0701 


5.5058 


3 


9637 


4 


7118 


5 


0477 


4 


3351 


3 


7336 


4 


3172 


4 


4860 


""Te 


PQQl 


5 


7440 


6.3018 


4 


6613 


5 


4025 


5 


7319 


5 


0177 


4 


4319 


5 


0103 


5 


1753 




PQQHHl 


4 


9231 


5.3418 


3 


8530 


4 


5817 


4 


9067 


4 


2084 


3 


6277 


4 


1964 


4 


3595 




PQQ2 


5 


6568 


6.2055 


4 


5875 


5 


3192 


5 


6446 


4 


9397 


4 


3610 


4 


9320 


5 


0951 




PQQHH2 


4 


9115 


5.3304 


3 


8459 


4 


5714 


4 


8951 


4 


1999 


3 


6218 


4 


1879 


4 


3503 


isoNd 


PQQl 


4 


1436 


4.5674 


3 


3937 


3 


9137 


4 


1420 


3 


6355 


3 


2316 


3 


6375 


3 


7514 




PQQHHl 


3 


1506 


3.4603 


2 


5501 


2 


9650 


3 


1478 


2 


7448 


2 


4208 


2 


7447 


2 


8359 




PQQ2 


4 


0499 


4.4632 


3 


3160 


3 


8249 


4 


0483 


3 


5526 


3 


1574 


3 


5546 


3 


6661 




PQQHH2 


3 


2415 


3.5638 


2 


6341 


3 


0545 


3 


2392 


2 


8305 


2 


5031 


2 


8311 


2 


9234 



B. Uncertainties in NTMEs 

To estimate the uncertainties associated with the 
NTMEs Af^"") for {P'p-)^^ decay calculated using the 
PHFB model, we evaluate their mean and the standard 
deviation, defined as 



(0<.) 



N 



(18) 



and 



1 



■ N 

E 

.4 = 1 



M. 



1/2 



(19) 



Recently, it has been shown by Simkovic et al. [8| that 
the phenomenological Jastrow correlations with Miller- 
Spenser parametrization is a major source of uncertainty. 
Therefore, it is more appropriate to consider SRC2 or 
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TABLE II: Maximum and minimum relative change in the NTME M^""^' (in %), for all nuclei included in table I, due to the 
use of a different average energy denominator ( second column), the inclusion of three different parameterizations of the SRC 
(SRCl, SRC2 and SRC3) with point nucleons (third to fifth column), the inclusion of finite size effect (F) (sixth column) and 
finite size effect plus SRC (F+SRCl, F+SRC2 and F+SRC3 in last three columns). In each row, the results employing one of 
the four different parameterizations of the effective two-body interaction are displayed. 



Parametrizatios 


A/2 




P+S 




F 




F+S 








SRCl 


SRC2 


SRC3 




SRCl 


SRC2 


SRC3 


PQQl 


7.9-10.2 


18.1-22.0 


5.5-7.3 


0.04-0.7 


12.3-13.9 


22.0-26.4 


12.2-15.0 


9.5-11.6 


PQQHHl 


7.5-9.8 


19.1-22.9 


5.9-7.7 


0.1-0.7 


12.9-15.0 


23.2-27.6 


12.9-15.7 


10.0-12.5 


PQQ2 


7.9-10.2 


18.1-22.1 


5.6-7.4 


0.04-0.7 


12.2-14.5 


22.0-26.6 


12.2-15.1 


9.5-11.7 


PQQHH2 


7.4-9.9 


18.7-23.1 


5.8-7.8 


0.1-0.8 


12.7-15.1 


22.8-27.8 


12.7-15.8 


9.8-12.2 



TABLE III: Average NTMEs AP and uncertainties AM* 
for the decay of ^^'^^Zr, ^^-I^OMo, ""Pd, i^s^i^OTe 

and ^^"Nd isotopes. Both bare and quenched values of qa are 
considered. 



ri3- 


9A 


Case I 


Case II 


emitters 










AM^"^' 


y^Zr 


1.254 
1.0 


4.2464 
4.6382 


0.3883 
0.4246 


4.4542 
4.8668 


0.2536 
0.2759 




1.254 
1.0 


3.1461 
3.4481 


0.2778 
0.3085 


3.3181 
3.6376 


0.1243 
0.1424 


9«Mo 


1.254 
1.0 


7.1294 
7.8398 


0.6013 
0.6826 


7.4656 
8.2099 


0.3635 
0.4358 


iooMo 


1.254 
1.0 


6.8749 
7.5660 


0.6855 
0.7744 


7.2163 
7.9419 


0.4977 
0.5769 


"«Pd 


1.254 
1.0 


7.8413 
8.6120 


0.8124 
0.9184 


8.2273 
9.0370 


0.6167 
0.7128 


i28Te 


1.254 
1.0 


4.0094 
4.4281 


0.4194 
0.4601 


4.2175 
4.6571 


0.3074 
0.3355 




1.254 
1.0 


4.4458 
4.9065 


0.5231 
0.5837 


4.6633 
5.1459 


0.4269 
0.4802 




1.254 
1.0 


3.1048 
3.4334 


0.4649 
0.5181 


3.2431 
3.5856 


0.4434 
0.4952 



SRC3 due to the Argonne VI 8 and CD-Bonn NN po- 
tentials, respectively. Based on these observations, we 
perform the statistical analysis of two cases. In case I, 
we calculate the average and variance of twelve NTMEs 
listed in the last three columns (F-l-S) of Table U with 
the bare and quenched values of axial vector coupling 
constant = 1.254 and gA = 1.0, respectively. The 
average and standard deviations of eight NTMES M^O") 
due to SRC2 and SRC3 are similarly calculated in the 

case IL The average NTMEs m''^'^'' and standard de- 
viations AM^*^"^ are presented in Table IIIII It is no- 
ticed that the exclusion of Miller-Spenser parametriza- 
tion reduces the uncertainty by about 55% in ^^Zr to 



4% in ^^°Nd isotope. In Table |TVl we present the av- 
erage NTMEs ' of case II along with the recently 
reported results in ISM by Caurier et al. QRPA as 
well as RQRPA by Simkovic et al. dl and IBM by Barea 
and lachello [31 . In spite of the fact that different model 
space, two-body interactions and SRC have been used in 
these models, the spread in the NTMEs turns out to be 
about a factor of 2.5. Further, we extract upper limits 
on the effective mass of light neutrinos {mu) from the 
largest observed limits on half-lives '^^ iP~ 
cay using the phase space factors of Boehm and Vogel 
p5| . It is observed that the extracted limits on (m^) 
for loOMo and ^^^Te nuclei are 0.48+°;°^ - 0.69+°;°^ and 

0-30to:o2 - 0-42lo:o4 eV, respectively. In the last column 
of Table HVl the predicted half-lives of (/3~/3~)g decay 
of 94,962r, 98'100mo, "Opd, i28,i30Te and i^^Nd isotopes 
are given for (m,y) — 50 meV. 



IV. CONCLUSIONS 

We have studied the (/3~/3~)oi/ decay of ^''^^^Zr, 
98,100^0, 104ru^ iiopd, I28,i30xe and ^^ONd isotopes in 
the light Majorana neutrino mass mechanism using a 
set of PHFB wave functions. The reliability of wave 
functions generated with PQQl and PQQHHl interac- 
tions has been tested in previous works by calculating the 
yrast spectra, reduced B{E2:0'^ — >• 2+) transition proba- 
bilities, static quadrupole moments Q(2+) and 5-factors 
(7(2+) of participating nuclei in P~)2i, decay as well 
as M21, and comparing them with the available experi- 
mental data [33, Hll . An overall agreement between the 
calculated and observed spectroscopic properties as well 
as M21, suggests that the PHFB wave functions generated 
by fixing Xpn or Xpp to reproduce the E2+ are reasonably 
reliable. 

In the present work, NTMEs M'-^"^ were calculated 
employing the PHFB model with four different parame- 
terizations of the pairing plus multipolar type of effective 
two body interaction and two (three) different parameter- 
izations of the short range correlations. It was found that 
the NTMEs Af^"") change by about 4-14(10-15)%. 

The mean and standard deviations were evaluated 
for the NTMEs M^"") calculated with dipole form fac- 



7 




FIG. 2: Radial dependence of C^°"\r) for the (/3"/3")(,^ decay of '"'Zr, "'°Mo, "''Pd, i^s.isorpg ^^^^^ isoj^^ isotopes. In this 
Fig., (a), (b), (c) and (d) correspond to P, P+SRCl, P+SRC2 and P+SRC3, respectively. Further, (e), (f), (g) and (h) are for 
F, F+SRCl, F+SRC2 and F+SRC3, respectively. 
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TABLE IV: Extracted limits on effective light Majorana neutrino mass (m^) and predicted half lives using average NTMEs 
m'""^ and uncertainties AM*""' for the {rDo, decay of ^^'^^Zr, ^^'^ooMo, "°Pd, ^^^-^^o^^ iso^^ isotopes. 



emitters 


OA 


-irT(Ov} 

M ' 


ISM 
[9] 


(R)QRPA 

m 


IBM 

m 


Goi 

(10-^V) 


Tlj^{ yr) Ref. 


{my) 


r{';2(y) 

< TTiu >= 50 meV 


^^^Zr 


1.254 
1.0 


4.45±0.25 
4.87±0.28 








0.1684 


1.9x10'" 


[46] 


&Alt',^ X 10^ 
9.23lg j^ X 10^ 


6.48t;;;^^xl02^ 


^«Zr 


1.254 
1.0 


3.32±0.12 
3.64±0.14 




1.43-2.12 




5.930 


1.0x10^' 


[46] 


2o.oo«:^^ 

28.70ll:J^ 


1.601°}? X 10^6 
3.29lgilxl02'^ 


9«Mo 


1.254 
1.0 


7.47±0.36 
8.21±0.44 








0.0018 


1.0x10" 


m 


imtlZ X io« 

2.32^0." X 10'^ 


1.06+o}Jxl02s 
2.16t°J5xl02« 


i°°Mo 


1.254 
1.0 


7.22±0.50 
7.94±0.58 




2.91-5.56 


3.732 


4.640 


4.6x10^^ 


m 


0.48«°* 
0-69t°°^ 


4.32toj|xl02^ 
8.83tiigXl02^ 


"«Pd 


1.254 
1.0 


8.23±0.62 
9.04±0.71 








1.422 


e.oxio''^ 


[49] 


6.72toi^ X 10^ 
9.63lo TO X 10^ 


1.09+°}«xl02« 
2.22t°;«xl02s 


i28Te 


1.254 
1.0 


4.22±0.31 
4.66±0.34 


2.26 


3.21-5.65 


4.517 


0.1849 


1.1x10^^ 


[50] 


8.50«JI 
12.10l°|t 


3.18+{jgxl02^ 
6.44lJ-:°lxl02^ 


"OTe 


1.254 
1.0 


4.66±0.43 
5.15±0.48 


2.04 


2.92-5.04 


4.059 


4.490 


3.0x10^" 


M 


sn+°"^ 

"•'JU_o.o2 
r, .0+0.04 


1.07+o-f,xl02« 
2.17t°;*^xl0^« 




1.254 
1.0 


3.24±0.44 
3.59±0.50 






2.321 


21.16 


1.8x10^^ 


[52] 


9 i;i:+0.40 

Q f^o+0.58 
O.00_o.44 


4.69li;ggxl0^^ 
9.49l^;?^xl0^'^ 



tor and with and without Miller- Spencer parametriza- 
tion of short range correlations, which were employed 
to estimate the )oi/ decay half-lives T^j2 for both 

QA = 1.254 and qa = 1-0. The largest standard devia- 
tion, interpreted as theoretical uncertainty, turns out to 
be of the order of 15% in the case of -'^^°Nd isotope. We 
have also extracted limits on the effective mass of light 
Majorana neutrinos (m^) from the available limits on ex- 
perimental half-lives T^j2 using average NTMEs M ^^"^ 

calculated in the PHFB model. In the case of ^^°Te iso- 
tope, one obtains the best limit on the effective neutrino 
mass {m^) < 0.30l[J;S]^ - 0.42+°;°^ eV from the observed 
limit on the half-lives T^J^ > 3.0 x 10^^ yr of (/3",9~)o^ 
decay fBT]- 

Note: Due to an error in one equation, the NTMEs 



Mp, Mqt, M^O"), Mph, MoTh andAPj^"' given m Ref. 
122] must be multiplied by a factor of 2. It implies that 
the limits on the effective light neutrino mass < > 
must be reduced by a factor of 2 whereas the limits on 
effective heavy neutrino mass < Mn > must be multiplied 
by a factor of 2. In both cases the limits are twice more 
stringent. 
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